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Abstract. For reproducing the anomalous — i.e., sub- or super-difFusive — behavior 
in some stochastic dynamical systems, the Generalized Langevin Equation (GLE) has 
gained considerable popularity in recent years. Motivated by the question whether or 
not a system with anomalous dynamics can have the GLE formulation, here I consider 
polymer physics, where sub-difFusive behavior is commonplace. I provide an exact 
derivation of the GLE for phantom Rouse polymers, and by identifying polymeric 
response to local strains, I argue the case for the GLE formulation for self-avoiding 
polymers and polymer translocation through a narrow pore in a membrane. The 
number of instances in polymer physics, where the anomalous dynamics corresponds 
to the GLE, thus seems to be fairly common. 
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Following the work of Einstein and Smoluchowski, a century ago Langevin proposed 
a stochastic dynamical description of Brownian motion [1]. The corresponding Langevin 
equation (LE) for a particle of mass m and velocity v{t) in a fluid of damping coefficient 
7 (meaning that the friction coefficient is rwy) is given by 



where the stochastic force f{t) satisfies (/(t)) = and the fluctuation-dissipation 
theorem (FDT) {f{t)f{t')) = 2m'jkBT6{t — t'). Here fc^ is the Boltzmann constant, 
T is the absolute temperature of the fluid and the angular brackets denote equilibrium 
ensemble averaging. Having found wide-ranging applications ^ the LE ([1]) — describing 
the Brownian motion in the so-called Rayleigh formulation — has long since been 
incorporated into the fundamentals of stochastic processes [31 H]. An exact result of 
the LE is that the dynamics of the Brownian particle is diffusive at long times, with the 
diffusion coefficient D = ksT / {ni'^) |5]. However, there is a diverse range of physical 
processes, for which the diffusive behavior at long times is not the norm; instead the 
mean-square-displacement (MSD) of the particle is anomalously fast or slow, increasing 
in time as for a ^ 1. For describing the dynamics in some of these systems, e.g., 
turbulent diffusion |6], disorder related excitations [7], ATP coupling to motor proteins 
[8], dipolar chains in a ferrofluid [9], ferrofluid patterns in a magnetic field p^, and 
traffic flows [11], a generalization of the LE in the Rayleigh formulation ([1]) — the 
so-called Generalized Langevin Equation (GLE) of the Mori-Lee form [12] — given by 



has gained considerable popularity in recent years. Here, the stochastic force g{t) 
satisfles {g{t)) = and the corresponding FDT {g{t)g{t')) = mkBTT[t — t'). The 
GLE reduces to the LE when T{t — t') = 2^6(t — t'). It is the non- instantaneous nature 
of r(t) in the GLE that leads to the anomalous dynamics: the result that if r(t) ~ 
for some a at long times, then the particle's MSD ~ t", has been derived not so long 
ago [13]. 

The characterization of the anomalous dynamics, given the characteristics of r{t) 
in the GLE, is a one-way street. From this perspective, it would also be worthwhile to 
know whether or not a system with anomalous dynamics can have the GLE formulation. 
While a generic answer to this question is not known to the best of my knowledge, I note 
that the complexity of the systems that exhibit anomalous dynamics typically presents 
a barrier for the answer. Motivated by this question, I consider in this Letter polymer 
physics, where anomalous dynamics is commonplace. I derive the GLE exactly for 
the motion of a tagged monomer for a phantom Rouse polymer. [The phantom Rouse 
equation follows the so-called Smoluchowski formulation [M], hence the corresponding 
GLE is not of the Mori-Lee form ([2]); see later]. Further, I identify polymeric response to 
local strains, and argue the existence of the GLE for (single) self-avoiding polymers and 
polymer translocation. Interestingly, despite the fact that any flrst course on polymer 
physics teaches that the dynamics of a tagged monomer of a polymer is anomalous until 



mv{t) 



mjv{t) + f{t), 



(1) 




(2) 
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the terminal time r (the relaxation time of the polymer), I have not seen, in published 
literature, the derivation of the GLE for polymer dynamics, and subsequently, the follow- 
up to the anomalous dynamics. This Letter thus provides the first indication that the 
number of instances in polymer physics, where the anomalous dynamics corresponds to 
the GLE, may indeed be fairly common. Whether or not this procedure for the GLE 
formulation can be extended to other systems (that exhibit anomalous dynamics) is also 
brought into consideration. 

The correct single polymer dynamics in a fluid was first presented by Zimm 
|15j . I refer to such polymers, for which the monomers interact with each other 
via hydrodynamic interactions, in this Letter, as Zimm polymers. Few years earlier 
than Zimm, Rouse [16] put forward a model for single polymer dynamics that 
neglect the hydrodynamic interactions between the monomers; although incorrect, the 
corresponding (Rouse) polymer dynamics resides at the heart of polymer physics [H] 
— and widely used till today — largely due to its simplicity. 

The GLE for phantom Rouse polymers. Consider a phantom Rouse polymer (of 
length N); i.e., a Rouse polymer that can intersect itself. It is described by the Rouse 
equation; in continuum representation it reads [11] 

where fn{t) is the location of the n-th monomer at time t, 7 is the damping coefficient 
of the surrounding fluid, and k is the spring constant for the springs connecting the 
consecutive monomers. The stochastic force /m(^) satisfies the conditions {fn(t)) = 
and {fmK(t)fn\{t')) = 2-fkBT6{m - n)6^xS(t - t'), for k, A = {x,y,z). Equation is 
supplemented by the "open" boundary conditions that the chain tension of the polymer 
at the free ends must vanish; i.e., {drn/dn)\n=o = {drn/dn)\n=N = 0. 

Since the Rouse equation is linear in rn{t), it can be solved to obtain all correlation 
functions using the mode expansion technique [H]. Two noteworthy results borne out 
of this exercise are: (a) the terminal relaxation time r = '-fN^/ln^k), and (b) the MSD 
of the middle monomer increases as t^^"^ until time r, and only after that time the middle 
monomer performs diffusive motion, with the diffusion coefficients scaling as It is 

this sub-diffusive motion of the middle monomer that I obtain from the GLE, which in 
turn I derive exactly from Eq. More precisely, I demonstrate that the chain tension 
(f){t) that the middle monomer experiences is obtained from the velocity of the middle 
monomer v{t) via the GLE 

${t) = - f dt' ^i{t-t')v{t')+g{t), (4) 



with g{t) = (j){t)^=o, so that {g{t)) = 0, and {g{t) ■ g{t')) = ^ksT ^i{t - t') = 
6kBTy/7nk{t - t'y^/^ exp[-{t - t')/T] is the FDT. The factor 3 in the expression for 
{g{t) - git')) stems from the fact that I consider polymers in three dimensions. From Eq. 
(jlj) I further arrive at the GLE 

v{t) = - / dt' a{t - t') (fit') + h{t), (5) 
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which encodes the anomalous dynamics of the middle monomer. Here {h{t)) = 0, 
{h{t) ■ h{t')) = SksT a{t — t') is the FDT, and fi{s)a{s) = 1 in the Laplace space. 

Two ingredients are necessary to derive Eq. (jl]). The first one of them is the 
dynamics for both halves of the polymer when the middle monomer held fixed at, say, 
n. With f'„(t) = fn{t) - n, I define 

— 

Y^'\t) = ^IJdn sin ^^i^E±}h ^^/^(,), (6) 

and Y(^\t) = -1 r dn sin ^li^Ellh _^^/^(,) (7) 

for p = 0, 1, . . ., for the right and the left half, such that 

rU./.(t) = 4^ffl(t)sin "^'^+^^" e(n) 



N 

V 



-45:y;<')Wsm^e^e(-„). (8) 

P 

Equation ([8]) shows that the open boundary conditions are satisfied, while the middle 
monomer remains fixed. The independent evolution of each half satisfies 

dY - - 

ip^ = -Kyvit)+fpit)^ (9) 

for % = [Yi'\Yi% 7p = 2iV7, = 2Ti''k{2p + if / N , (U) = and (/p«(t)/,A(t')) = 

^pq^KXlpkBTdit - 1'). 

The second ingredient is the equilibrium averages of Yp and F/. First, (Yi'\t)) = 

{Yi^\t)) = by isotropy. Secondly, by left-right symmetry {[Yi'\t)f) = {[Yi^\t)f); 
and they are obtained from the LE (Q as 

{[Yi'\t)r) = {[Yi'Ht)?) = 3NksT/[An'k{2p + in (10) 

With these ingredients, I now follow the dynamics of the middle monomer. Consider 
the case where the middle monomer of a polymer, stationary at t = 0~, moves by a 
distance Sro at t = 0, corresponding to v{t) = 6rQ6{t). Then, for Yp = [Yp''\ Yp'^'^], I have 

Yp{0+) = Yp{0-) - 5ro/[7r(2p + 1)] (11) 

Until the time the middle monomer moves again, it feels k{df'^j^j^i^/dn)n=Q+ and 
—k{df'^^jy^2/'^^)n=o- forces from the right and the left half of the polymer respectively. 
The total force (f){t) experienced by the middle monomer is their sum, and from Eq. ([8]) 
it is given by 

*) = 4^ E ^^^^ [^/^(^) + ■ ^^2) 

p 

To obtain the time evolution of (j){t) until the middle monomer moves again, we return 
to Eq. Qj and write 



1 /■* 

i;(t)=e-M/7py^(0+) + - / cit'e~'=-(*~*')/>/p(t') 



(13) 
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for Fp = [Fp^'^fp^"^], yielding 

f^->oE^-'"'^' ^^YP^m. with 



q{t)/k g{t)/k 

^^(t)=e-V/>[yW(o-) + i;(0(o-)] 

+ - fdt' e-^-(*~*')/^- [Z^'-) it') + it')] . 
Ip Jo 

It is now seen, by converting the sum to an integral, that 



(14) 



q{t) = -f 5ro J2 ^-''"'^ = -2^^-o \j ^ 6"*/^ (15) 
p 



It is also seen, using Eq. (fTOjl . that {git))^ = 0, with the FDT 



Here (• • .)o denotes the average over an equilibrated ensemble of polymers at t = 0~; 
for which the middle monomer of each polymer moves by 6ro at t = 0. 

The above procedure is trivially generalized to obtain the GLE one needs 
to consider an ensemble of polymers that moved by a distance 6ro at t = 0, 6ri at 
t = ti, and so on. For this ensemble, having recognized that v(t) = 'Ylii^'^i^i^ ~ ^i) 
[with to = 0], where the angular brackets include an average over an equilibrated 
ensemble of polymers at t = 0~, one arrives at Eq. (jl]). Further, the GLE ([5]) is 
obtained by first Laplace transforming Eq. (jl]), then expressing the velocity of the 
middle monomer in terms of the chain tension it experiences in the Laplace space, and 
finally inverting the Laplace transform to return to real time, resulting in the FDT 
{v{t) ■ t/(t'))^=o = ^ksTait - t') ~ (t - t')-3/2g-(t-t')/r Subsequently, the result 

that the MSD of the middle monomer increases ~ t^^^ till time r and ~ t thereafter is 
obtained by integrating of {v{t) ■ v{t'))^^Q twice in time. 

At this point I make the important observation for the GLE formulation of Eqs. 
(jlHSj) that if /i(t) ~ t""", then the MSD of the middle monomer has to increase ~ t". 

Polymeric response to local strains, and the GLE. From Eq. (jlj) it is clear that 
n{t) is the mean response of the polymer to a local strain — i.e., altered chain tension 

— * 

at the middle monomer — created by moving the middle monomer by a distance Sr 
at t = and fixing it at its new position Vt. The mean local strain then relaxes in 
time ~ t-i/^ i.e., (0(t) • 0(0)) ~ t-^/\ While the response of a polymer to a local 
strain depends on how the strain is created, the identification of fi{t) as the polymer's 
mean local strain relaxation response alone allows one to write down the GLE, as I 
show below. First, given the identification of fi{t) as the polymer's mean local strain 
relaxation response one can always write the stochastic Eq. (jll) with {g{t)) = 0, which 
holds by definition. Next, to obtain the FDT, consider Eq. (jlj) for an ensemble of 
polymers with v{t) = Wt and 4>{to) = go, a specific value. For such an ensemble 
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g[t) = (j){t), and since fi{t) is the polymer's mean local strain relaxation response, 
{4>{t) ■ 0(to)) = do f'it ~ to) for t > to- Extending this to the dynamics of a polymer 
in an equilibrium ensemble (where is also chosen from the equilibrium ensemble), 
one has {g(t) ■ g(to)) = {4>(t) ■ 0(t'))^=o = (0^(^))j/=o A^(^ ~ t') [for a phantom polymer 

{(j)\t))^=0 = SkBT]. 

Using the identification of fi{t) in Eq. (j4]) as the mean polymeric response to 
local strain leading to the GLE, I now argue the existence of the GLE for self-avoiding 
polymers and for polymer translocation. 

The GLE for self-avoiding polymers. The monomers of a self-avoiding polymer 
interact over a long-range, which prohibits one from writing down an exact equation for 
the velocities of the monomers in terms of the forces they experience. However, quite a 
few properties of self-avoiding polymers are well-known: two of them we need here for 
a polymer of length are: (i) the terminal time r scales ~ A^^^^'^ for a Rouse polymer, 
and as ~ A^^'^ for a Zimm polymer [13] ; and (ii) the entropic spring constant of a polymer 
scales as N~^'^ [19]. Here u is the Flory exponent, in three dimensions u ~ 0.588, and 
in two dimensions u = 0.75. Imagine that one moves the middle monomer of a self- 
avoiding polymer by a small distance 6r a.t t = and holds it at its new position, 
corresponding to v{t) = 6r6{t). Following (i), at time t, counting away from the middle 
monomer, all the monomers within a backbone distance Ut ~ ti/(i+2i') f^j. Rouse, and 
~ ^1/(3!^) for a Zimm polymer equilibrate to the new position of the middle monomer. 
However, since the rest (A^ — Ut) monomers are not equilibrated to the new position of 
the middle monomer at time t, these nt monomers are stretched by a distance 6r. With 
the entropic spring constant of these Ut equilibrated monomers scaling ~ 'n,^^'' [following 
(ii)], the mean force the middle monomer will experience at its new position is given by 
(f{t) ~ ut^^'i-d^r) ~ t-2'^/(i+2'^)(-5r) for a Rouse, and ${t) ~ n^'^'^i-d^r) ~ t-^/^{-6^r) 
for a Zimm polymer [force = (spring constant) x (stretching distance)]. This power-law 
behavior lasts only till the terminal time r. [The time behavior of Eq. ( |T5|) is recovered 
from this line of argument upon simply replacing u hj 1/2 corresponding to a phantom 
Rouse polymer.] Such behavior of yu(t), in light of the above paragraphs, implies that 
the motion of the middle monomer of the Rouse and the Zimm polymers is indeed 
described by the GLE, reproducing the well-known results that the MSD of the middle 
monomer increases respectively ~ t2!./(i+2!.) ^^^^ ^ ^2/3^ ^^^^ ^ ^ thereafter. The GLE 
for a self-avoiding Rouse polymer has recently been confirmed numerically [18] . 

The GLE for polymer translocation. Polymer translocation is the process where a 
polymer passes through a narrow pore in a membrane. Of interest here is the so-called 
unbiased (i.e., in the absence of any force or field) translocation: the polymer passes 
through the pore purely due to thermal fluctuations, and the dynamics is anomalous 
|20j . A translocating polymer essentially consists of two sub-polymers — one on each 
side of the membrane — exchanging monomers through the pore. When a monomer 
translocates, the polymer's chain tension at the pore changes: it increases on the side 
of the membrane which the monomer comes from, and decreases on the other. The 
relevant polymeric response therefore, is to a (local) strain due to extra monomers 
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injection into a tethered polymer at the tether point. Consider the case where n extra 
monomers are injected into a tethered polymer at the tether point at t = 0. For phantom 
Rouse polymers the mean response to such a strain is given by fi{t) ~ t-^e"*/^, with 
r ~ N"^ [21J. For self-avoiding polymers is obtained as follows. Following (i), at 
time t, counting away from the pore, all the monomers within a backbone distance 
Uf. ~ ^1/(1+2*^) for a Rouse, and ~ ^1/(3*^) for a Zimm polymer, equilibrate to the 
injected monomers. The real space extent of Ut monomers is r{nt) ~ n^, but since 
the rest (A^ — Ut) monomers are not equilibrated to the injected monomers at time 
t, there are {uf + n) monomers squeezed in a space that extends only to r{nt). The 
corresponding compressive force [force = (spring constant) x (stretching distance)] 
from these (rif + n) monomers, felt at the pore, and hence /i(t), is the given by 
~ n^'^'^[6r{nt)] ~ n^'^" n[dr{nt)/dnt] = unn^ ^^'^'^\ which scales ~ j(-(i+'^)/(i+2i') for a 
Rouse and ~ t-(i+^)/(3^) for a Zimm polymer. (Once again, this behavior lasts only 
till the terminal time r.) This implies that polymer translocation is described by the 
GLE as well, resulting in the scaling of the MSD ~ t(i+'^)/{i+2i') for self-avoiding Rouse 
and ~ for self-avoiding Zimm polymers up to time r and ~ t thereafter. 

Consequently, the pore-blockade time scales ~ A^^ for a phantom Rouse [21], ~ N'^'^'^ 
for self-avoiding Rouse [221 [23], ^ind ~ for a self-avoiding Zimm polymer [22j . 

In summary, in view of the recent popularity of the GLE to reproduce the anomalous 
dynamics in some stochastic dynamical systems, in this Letter I have concerned myself 
with the question whether a system with anomalous dynamics can lead to the GLE 
formulation, and have considered polymer physics, where sub-diffusive behavior is 
commonplace. I have provided an exact derivation of the GLE for phantom Rouse 
polymers, and have argued the case for the (necessary) existence of the GLE for self- 
avoiding polymers and for polymer translocation. The key to show the existence of 
the GLE for these cases is to relate the (monomeric) velocity to the power-law mean 
(polymeric) response behavior of to local strains. This implies that there exists the GLE 
formulation for any system (that exhibits anomalous dynamics) for which the velocity 
gives rise to local strains, and the local strain relaxes in a power-law fashion. Moreover, 
all the cases considered here concern single polymer dynamics; whether the principle 
holds for many-polymer systems, e.g., for subdiffusive behavior in polymer melts, or for 
super diffusive behavior in "living polymers" [25] remains to be investigated. 
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